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ABSTRACT 

We analyze the long-term (rest-frame 3-30 yr) X-ray variability of eleven broad absorption line 
(BAL) quasars, mainly to constrain the variation properties of the X-ray absorbing shielding gas 
that is thought to play a critical role in BAL wind launching. Our BAL quasar sample has coverage 
with multiple X-ray observatories including Chandra, XMM-Newton, BeppoSAX , ASCA, ROSAT, 
and Einstein; 3-11 observations are available for each source. For seven of the eleven sources we 
have obtained and analyzed new Chandra observations suitable for searching for any strong X-ray 
variability. We find highly significant X-ray variability in three sources (PG lOOH-054, PG 1004-1-130, 
and PG 2112-1-059). The maximum observed amplitude of the 2-8 keV variability is a factor of 
3.8 ± 1.3, 1.5 ± 0.2, and 9.9 ± 2.3 for PG 1001-h054, PG 1004-f 130, and PG 2112-^059, respectively, 
and these sources show detectable variability on rest-frame timescales down to 5.8, 1.4, and 0.5 yr. 
For PG 1004-1-130 and PG 2112-f 059 we also find significant X-ray spectral variability associated with 
the flux variability. Considering our sample as a whole, we do not find that BAL quasars exhibit 
exceptional long-term X-ray variability when compared to the quasar population in general. We do 
not find evidence for common strong changes in the shielding gas owing to physical rearrangement or 
accretion-disk rotation, although some changes are found; this has implications for modeling observed 
ultraviolet BAL variability. Finally, we report for the first time an X-ray detection of the highly 
polarized and well-studied BAL quasar IRAS 14026-1-4341 in its new Chandra observation. 
Subject headings: cosmology: observations — X-rays: galaxies — galaxies: active — quasars: absorp- 
tion lines 



1. INTRODUCTION 

Fast outflows are common in active galactic nuclei 
(AGNs) ov er a range of more th an 10,000 in AGN lumi- 
nosit y (e.g.. lLaor fc Bra ndt 2002: lGangulv fc BrothertonI 
12008). These outflows are a significant component of the 
nuclear environment, and their ubiquity suggests that 
mass ejection in a wind is fundamentally linked to mass 
accretion onto a supermassive black hole (SMBH). In- 
deed, outflows may be key in allowing the accreting mat- 
ter to expel angular momentum. Observationally, while 
these outflows are most commonly studied via the ab- 
sorption features they imprint upon ultraviolet (UV) and 
X-ray spectra, they are also thought to shape emission- 
line and continuum prope rties (e.g., iMurrav fc Chian j 
Il997t [Richards et al.|[201l[) . Studies of the outflows from 
luminous AGNs, such as those seen in Broad Absorp- 
tion Line (BAL) quasars, have taken on further im- 
portance in recent years as they are one of the possi- 
ble agents of feedback in typical massive galaxies (e.g., 
iDi Matteo et al.ll2005l ). Outflows from luminous AGNs 
may expel enough gas from the SMBH/galaxy system 
to quench both SMBH growth and star formation, per- 
haps leading to the obser ved relations betwee n SMBH 
and bulge properties (e.g.. iGiiltekin et al]l2009( ). 

One frequently used and generally well-supported 
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model proposes that the UV absorption lines of BAL 
quasars originate in a line-driven wind that is launched 
from the accretion disk at 10^^~^'' cm from the SMBH 
(e.g., IMurrav et all 119951 : iProga et all l2000l ). In this 
model, the weak X-ray emission seen from BAL quasars 
(e.g., Gallagher et al. 2006; Gibson et al. 2009) is due 
to absorption in stalled, highly ionized "shielding gas" 
at smaller radii that protects the wind from nuclear ex- 
treme UV and soft X-ray radiation; without such an ab- 
sorbing layer, these energetic photons would so highly 
ionize the potential wind material that it could not be 
effectively accelerated via line pressure. Supporting ob- 
servational evidence for the physical importance of the 
shielding gas comes in the form of correlations between 
the X-ray weakness of BAL quasars and the strength and 
velocity of their Uy C IV BALs (e.g.. iGallagher et all 



[200l[Gibson et al.ll200llWu et al.ll2CiM 

The nature of the shielding gas in BAL quasars 
remains poorly constrained observationally. Heavy 
X-ray absorption is directly seen in many of the BAL 
quasars with X.-ray spectra of reasonable q uality (e.g., 
'GaUa gher et"aI1 [2OOII . l2(¥)l iFan et al.ll2009h . Measured 
absorption column densities span a broad range of typ- 
ically Nn w 10^^-10^'* cm~^, although most individual 
iVjj measurements have significant uncertainties owing to 
limited photon statistics and absorption complexity (e.g., 
ionization and partial-covering effects). The kinematics 
of the X-ray absorber in BAL quasars is also poorly con- 
strained. The X-ray absorber is sometimes envisioned 
to be a stalled "sacrificial wind" that prot ects the UV- 
absorbing wind launched at larger radii ('e.g.. lProga et al.l 
I2000D . although in one BAL quasar apparent evidence 
for a rapidly o utflowing X-ray abso rber has bee n found 
fe.g.. IChartaset all 120091: ISaez fc Chart as! 1201 Ih . Mag- 
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netocentrifugal effects, as well as radiative ones, may play 
an important role in setting the kinematics of the X-ray 
absorbing shielding gas fe.g.. lEverettll2005l ). 

Remarkable X-ray absorption varia bility has been 
seen i n a few BA L qua s ars (e . g.. [Gallagher et alJ 
2001 iSchartel et alJ [20051 [20071: iMiller et alJ 120061 : 



Chart as et al. II2009D, typically on multi-year timescales. 



indicating that the X-ray absorber is indeed a dynam- 
ical structure. This variability could be due simply to 
accretion-disk rotation moving an azimuthally asymmet- 
ric X-ray absorber through our line of sight; at the ex- 
pected radius of the shielding gas (< 10^^~^^ cm), sub- 
stantial fractional disk rotation is expected on multi-year 
timescales. The results on X-ray absorption variability in 
BAL quasars are broadly reminiscent of those found for 
lower luminosity Seyfert 2 galaxies, where X-ray absorb- 
ing column densities are found to vary commonly (e.g., 
iRisaliti et al.l l2002l I2009D . In addition to variations in 
X-ray absorber column density, some BAL quasars with 
multiple observations have also shown notable continuum 
luminosity variability. For example, the normalization of 
the power-law continuum in PG 2112+0 59 decreased by 
a fac tor ss 3.5 between 1999 and 2002 fjGallagher et al.l 
|2004|) . and continuu m variability of this m agnitude con- 
tinued through 2007 ([Schartel et al.ll2010[ ). 

Current studies of the X-ray absorption variability of 
BAL quasars suffer from small sample sizes. Presently 
only a handful of BAL quasars have sensitive coverage 
in multiple X-ray observations, and thus the general fre- 
quency of BAL quasar X-ray variability remains poorly 
characterized. In order to improve this situation, we 
targeted seven optically bright BAL quasars as part of 
Chandra Cycle 11 that already had sensitive archival 
X-ray observations with Chandra, XMM-Newton, or ear- 
lier missions. We obtained relatively short observations 
(5-7 ks) suitable for searching for any strong X-ray vari- 
ability. Our main goal was to execute an economical 
multi-year X-ray variability search for as many objects 
as possible, so that the basic X-ray variability properties 
of BAL quasars could be characterized as a class. Simi- 
lar ongoing efforts are also now characterizing the niulti- 
year UV absorption - line v a riability of BAL quasar s (e.g. , 
2?)^ [2OIOI : iCapelluDO et~aIll20Tll [2OII 



Gibson et al.l 



Filiz Ak et al 



I2OI2I ) , and X-ray constraints upon shield- 
ing gas variations will inform these efforts. In addition 
to presenting our new Chandra Cycle 11 observations be- 
low, we also utilize the existing archival X-ray observa- 
tions of BAL quasars capable of constraining long-term 
X-ray variability; this maximizes our sample size. 

The layout of this paper is as follows: in we de- 
scribe the sample and the X-ray data reduction; in ^ 
we provide a description of the spectral analysis; in ^ 
we provide our variability results and discussion focusing 
on long-term X-ray variability; and in ^ we summarize 
our main findings. Throughout this paper, unless stated 
otherwise, we use cgs units, the errors listed are at the 
la level, and we adopt a flat A-dominated universe with 
Ho = 70 km s-i Mpc"\ ilA = 0.7, and = 0.3. 

2. OBSERVATIONS AND DATA REDUCTION 

Our sample consists of eleven BAL quasars that have 
been observed in the past by sensitive X-ray missions 
including Chandra, XMM-Newton, BeppoSAX, ASCA, 
ROSAT, and Einstein. Seven of these BAL quasars 



have been re-observed with new short-exposure (5-7 ks) 
Chandra observations that are analyzed in this paper. 
In Table [T] we present the Galactic column densities, 
radio- loudness parameters (Rl),^ and some basic optical 
and radio properties for our sample. The redshifts of our 
sample quasars range from 0.04—2.88. The only sources 
that are not radio-quiet quasars (RQQs; Rl < 10) are 
PG 1004-^130 and Mrk 231, radio-loud quasars (RLQs) 
with Rl > 10. Notice, however, that in the case of 
Mrk 231 there is likely substantial reddening {Ay ~ 2; 
e.g.. lBoksenberg et al.il977i) : as a consequence the intrin- 
sic 4400 A flux density would be approximately an order 
of magnitude higher. Based on this argument Mrk 231 
would not be considered a bona- fide RLQ. 

In Table [T] we also present the number of X-ray obser- 
vations of each source and the rest-frame time span cov- 
erage. The distributions of the number of X-ray observa- 
tions with the respective rest-frame time span coverage 
are presented in Figure [TJ There are 3—11 X-ray ob- 
servations with rest-frame time span ranging from 3—30 
yr per source. The observation logs of the sources an- 
alyzed in this work, which include observation dates, 
observed count rates with Icr error bars, and exposure 
times, are presented in Table [5] for the Chandra, Bep- 
poSAX, ROSAT, and Einstein observations; Table [3] 
for the XMM-Newton observations; and Table [H for the 
ASCA observations. The errors on the source counts 
were computed by propagating the asymmetric errors on 
the total and background co unts using the approach de- 
scribed in §5 of lBarl ow (2004). The total and back ground 
count errors were obtained from Tables 1 and 2 of iGehrelsl 
(1198^ . 

2.1. Chandra observations 

The Chandra observations (new and archival) were an- 
alyzed using the standard software CIAO version 4.2 pro- 
vided by the Chandra X-ray Center (CXC). Standard 
CXC threads were employed to screen the data for sta- 
tus, grade {ASCA grade 0, 2, 3, 4, and 6 events), and 
time intervals of acceptable aspect solution and back- 
ground levels. For most of the observations the source 
and background spectra were extracted (using the CIAO 
task dmextract) from a circular region with an aper- 
ture radius of 4" ^ and an annular source-free region with 
an inner radius of 6" and an outer radius of 24", re- 
spectively. For the cases of PG 1004+130 and Mrk 231 
the source, annulus inner, and annulus outer radii are 
(16", 20", 30") and (24", 30", 50"), respectively. For 
these sources we have selected larger extraction regions 
given that they possess extended X-ray emission. The 
extended emissi on of PG 1004+1 30 a nd Mrk 231 has 
been s tudied bv lMiller etall (|2006D and [Gallagher et al.l 
(|2002| ). respectively. We selected the radius of extraction 
for these sources using a fixed background region (annu- 
lar region of inner radius 40" and outer radius 60") and 
estimating the total counts as a function of aperture ra- 
dius. The aperture radius selected corresponds to that 
where the total counts reach an approximately constant 

® The radio-loudness parameter is defined as the ratio between 
the flux density at 5 G Hz and the flux de nsity at 4400 A in the 
rest frame of the source l IKellermann et aI.| [T989) . 

^ The 1.5 keV encircled energy at 4" from the center of an on-axis 
point source is ~ 99% for ACIS-S . 



3 



TABLE 1 

Optical/Radio Properties, Galactic Column Densities, and X-ray Coverage of Surveyed BAL Quasars 



OBJECT NAME 


z 


"2000.0'' 


"52000.0^ 


A^HGal'' 


ABasoo'^ 


log Ri^'^ 


BAL type° 






IRAS 07598+6508 


0.148 


08 


04 


30.51 


+64 


59 


52.7 


4.3 


14.5 


0.80 


Lo 


5 


16.8 


PC 0946+301 


1.221 


09 


49 


41.10 


+29 


55 


19.2 


1.7 


16.2 


-1.71 


Hi 


3 


7.3 


PC 1001+054 


0.160 


10 


04 


20.14 


+05 


13 


00.5 


2.4 


16.6 


-0.12 


Hi 


5 


26.0 


PC 1004+130 


0.241 


10 


07 


26.11 


+12 


48 


56.1 


3.7 


15.3 


2.18 


Hi 


3 


19.9 


Q 1246-057 


2.247 


12 


49 


13.86 


-05 


59 


19.1 


2.1 


16.3 


<-0.14 


Hi 


8 


7.7 


Mrk 231 


0.042 


12 


56 


14.23 


+56 


52 


25.2 


1.3 


15.3 


1.41 


Lo 


11 


28.4 


IRAS 14026+4341 


0.323 


14 


04 


38.80 


+43 


27 


07.4 


1.2 


16.2 


-0.27 


Lo 


4 


13.6 


CSO 755 


2.882 


15 


25 


53.89 


+51 


36 


49.2 


1.6 


16.2 


<-0.19 


Hi 


5 


2.9 


SBS 1542+541 


2.370 


15 


43 


59.45 


+53 


59 


03.4 


1.3 


16.8 


-0.85 


Hi 


4 


2.5 


PG 1700+518 


0.292 


17 


01 


24.80 


+51 


49 


20.0 


2.6 


15.3 


0.38 


Lo 


6 


15.0 


PG 2112+059 


0.466 


21 


14 


52.57 


+06 


07 


42.5 


6.6 


15.0 


-0.56 


Hi 


9 


10.9 



^ Optical positions in J2000.0 equatorial coordinates. 

Galactic absorption column densities in units of 10 cm' ffrom lDickev fc Lockmanll99(J) . obtained using the optical coordinates 
of the sources and the HEASARC A'^h tool at http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 

The monochromatic AB magnitude at rest-frame wavelength 2500 A. These were computed from the flux densities obtained from 
the closest optical magnitude found in the NASA Extragalactic Database (NED; http://ned.ipac.caltech.edu) and extrapolating a 
power- law with a = -0.5 (e.g. . IVanden Berk et al.ll200ll) . Previous to the extrapolation the flux densities obtained through NED 
have been corrected for Galactic reddening. 

Logarithm of the radio-loudness parameter {Rl — /^GHz/Z^^qq^)- The flux densities at rest-frame 4400 A are obtained using 
the Galactic reddening corrected flux densities obtained from the closest optical magnitudes and extrapolating a power-law with 
a — —0.5 fe.g.^ lVanden Berk et al.||200ll) . The flux densities at rest-frame 5 GHz have been obtained using either the flux densities 
at 1.4 GHz or the flux densities at 5 GHz. When both the 1.4 GHz and 5 GHz flux densities are available, we flnd the flux 
density at rest-frame 5 GHz through a power-law interpolation between the two values. For all the sources with the exception of 
SBS 1542+541 the flux densities at 1.4 GHz were obtained from the NVSS jGondon ct al. 19 91) the FIRST survey jWhite et al.| 
fT997h . Additionally for PG sources the flux densities at 5 GHz were obtained from Kellermann et al ] U98gD . For SBS 1542-^541 
the flux densities at 1.4 GHz and 5 GHz were obtained from lPetric et al.l II2006D . 

^ Hi — ultraviolet spectra show high-ionization BALs only; Lo — low-ionization (Mg II and/or Al III) BALs present. 
^ Number of X-ray observations (.^X-ray) with their rest-frame time span coverage in years (Aircst). 



value. The use of larger source regions for the Chandra 
observatfons of PG 1004-1-130 and Mrk 231 is necessary 
in order to make fair flux comparisons with the other X- 
ray missions that have more limited angular resolution. 
The numbers of counts in the source regions at ener- 
gies of 0.5-8 keV are presented in Table [2j The response 
files used to analyze the source spectra were created with 
the CIAO tasks MKACISRMF for the redistribution matrix 
files (RMFs) and mkarf for the ancillary response files 
(ARFs). 

2.2. XMM-Newton observations 

The XMM-Newton data were analyzed with the Sci- 
ence Analysis Software (SAS) version 10.0.0 provided 
by the XMM-Newton Science Operations Centre (SOC). 
The events files were obtained by including the events 
with flag = 0, pattern < 12 (pattern < 4), and 200 < 
PI < 12,000 (150 < PI < 15,000) for the MOS (pn) 
detectors. The events files were also temporally filtered 
in order to remove periods of flaring activity. To filter 
the events files for periods of flaring activity, we created 
light curves using single events with energies greater than 
10 keV. Using these light curves we find and remove 
from the events files the time periods where the count 
rate is above 0.5 count s~^ for the MOS cameras and 
1 count s~^ for the pn camera. For the three EPIC de- 
tectors (MOSl, M0S2, and pn), the source counts for 
each quasar were extracted from a circular region with 
an aperture radius of either 18" or 30".^ The default 
aperture radius of the source extraction region was 30"; 

® The 1.5 keV encircled energy at 18" (30") from the center of an 
on-axis point source is 75% (~ 85%) for the EPIC cameras. The 
aperture corrections for the EPIC cameras are performed through 
the use of the ARFGBN routine. 



however, when the background counts were > 20% of the 
total counts inside this region an aperture radius of 18" 
was used. Under this criterion most observations were 
reduced using 18" apertures. The cases that were re- 
duced with 30" apertures were the 2001 December 8th 
observation of CSO 755, the Mrk 231 observation, the 
PG 1001+054 observation, and the 2003 May 14 obser- 
vation of PG 2112-1-059. The background spectra for 
most of the sources were extracted from annular source- 
free regions with 50" inner radii and 100" outer radii. 
In cases where this region covered other X-ray sources 
or detector gaps, the background was extracted from a 
source-free circular region with radius 80" in the vicin- 
ity of the source. The numbers of counts in the source 
regions for the XMM-Newton cameras are presented in 
Table |3l These counts are selected at energies of 0.5- 
10 keV for the MOS cameras and 0.3-10 keV for the pn 
camera. The response files used to analyze the source 
spectra were created with the SAS tasks rmfgen for the 
RMFs and ARFGEN for the ARFs. 

2.3. ROSAT observations 

The ROSAT PSPC observations were analyzed follow- 
ing the ROSAT xselect guide. ^ Using the prescription 
given we applied all the standard corrections (gain varia- 
tions, gain saturation, analog-to-digital nonlinearity, and 
dead time) to obtain cleaned events which were used for 
the subtraction of the source and background spectra. 
The source spectra were extracted from circular regions 
with aperture radii of 50". The background spectra for 
most of the sources were extracted from annular source- 
free regions with 100" inner radii and 200" outer radii; 

http:/ /hcasarc. nasa.gov/docs/rosat/ros_xselect_guide 



TABLE 2 

Log of Chandra, BeppoSAX , ROSAT, and Einstein Observations 







EXP. TIME^ 




RATE 




OBJECT NAME 


OBS. DATE 


(ks) 


COUNTS^'*^ 


(10-3 s-1) 


Ref.'^ 




New Chandn 


1 observations (0.5-8 keV) 






IRAS 07598+6508 


2010 Jun 18 


6.7 






1 


PG 0946+301 


2010 Jan 11 


6.5 


57+9 


8.8ii:l 


1 


PG 1001+054 


2010 Jan 11 


4.9 


13li 




1 


Mrk 231 
IRAS 14026+4341 


2010 Jul 11 
2010 Jul 28 


4.8 
6.7 


247l« 


51.6+3'? 


1 
1 


CSO 755 
PG 1700+518 


2010 Mar 21 
2010 Jun 26 


4.9 
6.7 




22.4+2'5 


1 

1 



Archival Chandra observations (0.5-8 keV) 



IRAS 07598+6508 


2000 Mar 21 


1.3 


9+4 




2 


PG 1004+130 


2005 Jan 5 


41.1 


1905+1^ 


46.4+1.1 


3 


Q 1246-057 


2000 Feb 8 


5.5 




7g+1.4 
'•°-1.2 


2 


Mrk 231 


2000 Oct 19 


39.3 


2059^52 


52.5+1.3 


4 


Mrk 231 


2003 Feb 3 


39.7 


1947t*^ 


49.1+1.2 


5 


Mrk 231 


2003 Feb 11 


38.6 


1900+^^ 


49.2+1.2 


5 


Mrk 231 


2003 Feb 20 


36.0 




48.7+1.3 


5 


SBS 1542+541 


2000 Mar 22 


4.5 






2 


PG 2112+059 


2002 Sop 1 


56.9 


836129 


14.7+0.5 


6 


Archival BeppoSAX observations (1.8—10 keV) 


Mrk 231 


2001 Dec 29 


144.0 


818+36 


5.7+0.3 


7 


CSO 755 


1999 Feb 02 


35.2 


53+18 


1.5+0.5 


8 




Archival ROSAT observations 


(0.5-2 keV) 






IRAS 07598+6508 


1991 Mar 17 


7.9 




2.2+0.8 


9 


PG 0946+301 


1993 Nov 14 


11.8 


<7 


<0.6 


1 


PG 1001+054 


1992 May 17 


8.6 


Q+5 


1.1+0.5 


10 


Q 1246-057 


1991 Dec 21 


12.8 


20±8 


1.6+0.6 


1 


Q 1246-057 


1992 Jul 17 


2.8 


<4 


<1.5 


11 


Q 1246-057 


1993 Jan 5 


9.1 


14+^ 


^■"-0.5 


11 


Q 1246-057 


1993 Jun 22 


43.6 


5911^2 


1.4+0.3 


11 


Mrk 231 


1991 Jun 7 


23.2 


2241^^ 


9.7+0.8 


12 


IRAS 14026+4341 


1992 Jul 18 


6.2 


<5 


<0.8 


1 


SBS 1542+541 


1993 Ago 14 


5.5 






11 


PG 1700+518 


1991 Feb 9 


7.8 


<7 


<0.9 


11 


PG 2112+059 


1991 Nov 30 


11.7 






13 


Archival Einstein 


observations (0.16—3.5 keV) 




PG 1001+054 


1979 Dec 02 


1.6 


<15 


<9.4 


14 


PG 1004+130 


1979 May 21 


7.1 


<37 


<4.9 


14, 15 


Q 1246-057 


1979 Jul 20 


3.5 


<24 


<7.0 


14 


Mrk 231 


1980 Dec 23 


2.1 


<19 


< 8.1 


14 



^ The counts correspond to the net photon counts {background subtracted). The exposure 
times and photon counts arc obtained after screening the data. 

^ For details about the aperture sizes used see ^ The errors on the source counts were 
computed by propagating the asymmetric errors on the total and background counts using 
the approach of IBa rlow (2004). The total and background count errors were obtained from 
Tables 1 and 2 of lCehrcls 1 19861) . 

References: (1 ) This work; (2) Green et al. (2001); (3) Miller et al. (2006); 
(4) IGallagher et all | |2002|); (5) IGalla glicr et al. (2005); (6) Gallagher et al. (2004); 
(7)'Braito ct al. (2004); (8)IBrandt et all 11999): (9) Lawrence et al. (1997); (10) Laor et at] 
il997); (11) Green & Mathurl JToOl): fl2),R ,iggpoulou et al.. (.1996. '!: (13) . Wang et al. 1199a) : 
(14) IWilkes et all 1 119941) ; (15) IElvis fc Fabbiang (119841) . 
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TABLE 3 

Log of XMM-Newton Archival Observations 







EXPOSURE TIME^ 


MOSl'' 


MOS2'' 












MOSl MOS2 pn 


COUNTS RATEXlO'^ 


COUNTS RATE X 10'^ 


COUNTS 


RATE X 10'' 




OBJECT NAME 


OBS. DATE 


(ks) 


(0.5-10 keV) 


(0.5-10 keV) 


(0.3- 


-10 koV) 


REF.'^ 



IRAS 07598+6508 


2001 Oct 24 


18.6 


18.5 


14.2 


8411? 


4.5+0.6 


70+11 
"J- 10 


0+0.6 

-*-»-0.5 


268j;?9 


18.9i};* 


1 


PG 1001+054 


2003 May 04 


11.3 


11.3 


8.8 




6.9+0.9 




^+0.9 


322+20 


36.4j:^;3 


2 


PG 1004+130 


2003 May 04 


21.5 


21.5 


18.5 


567^27 


26.4+1.2 


622^2? 


28.9ti;2 


1777^1^ 


98.3i|« 


3 


Q 1246-057 


2001 Jul 11 


40.0 


40.0 


33.5 


126ti^ 


3.2+0.4 


127ti| 


3.2+0.4 


335+28 


10.0+0.8 


4 


Q 1246-057 


2004 Jul 15 


46.4 


47.0 


34.8 


141+16 


■^•'-'-0.3 


123ti« 


2.6+0.3 


351+25 


10.1+0.7 


5 


Mrk 231 


2001 Jun 07 


21.4 


21.4 


17.2 


774^2!) 


36.1+1.4 


739^29 


34.5+1.4 


2165t^J 


125.8+2.9 


6 


IRAS 14026+4341 


2005 Nov 26 


8.7 


7.9 


0.0 


<10 


<1.2 


14+6 


1 a+0.8 






7 


IRAS 14026+4341 


2006 Jan 16 


6.9 


6.2 


0.0 


12+6 


1 7+0.9 
^•'-0.8 


<8 


<1.3 






7 


CSO 755 


2001 Jul 30 


2.2 


2.2 


0.0 




13.11^;^ 


4lt? 








8 


CSO 755 


2001 Dec 08 


30.1 


30.3 


24.6 


610^2? 


20.3+0.9 


578^2^ 


19.1+0.9 


1820+^^ 


73.8+1.9 


8 


CSO 755 


2001 Dec 13 


4.5 


3.7 


2.0 


74+10 


16.3+21 


57+1 


15.2+2 2 


120+}* 


6i.ot^:? 


8 


SBS 1542+541 


2002 Feb 03 


8.4 


8.5 


2.4 


99+i? 


^^^•'-1.3 


74+11 
'^-10 




98ti? 


40.41^-0 


4 


SBS 1542+541 


2002 Feb 06 


23.8 


23.6 


12.7 


26lti^ 


11.0+0.8 


268^^8 


11.3+0.8 


48912^ 


38.5+1.9 


4 


PG 1700+518 


2009 Dec 19 


7.0 


6.9 


4.5 


<7 


<0.9 


<8 


<1.2 


30+9 


°-'-1.9 


9 


PG 1700+518 


2009 Dec 31 


14.1 


12.8 


7.9 


<12 


<0.9 


14+7 


1 1+0.6 

-'-•-'--0.5 


53ti? 


6 7+^-5 


9 


PG 1700+518 


2010 Jan 02 


19.3 


18.8 


15.4 


32t^ 


1 7+0.5 

-OA 


20+8 


1.1+0.4 


87+14 


5.6+0.9 


9 


PG 2112+059 


2003 May 14 


9.3 


9.7 


6.9 








35.3+2.1 


900+1* 


iso.ott^ 


10 


PG 2112+059 


2005 Nov 20 


72.9 


73.1 


64.8 


245+19 


3.4+0.3 


235ti^ 


r, + 0.3 


78511* 


12.1+0.5 


10 


PG 2112+059 


2007 May 03 


23.7 


23.2 


19.3 




5.6+0.6 


124+13 


c; O+0.6 
3-^-0.5 


336+23 


17.4+1.2 


10 


PG 2112+059 


2007 May 19 


59.3 


74.2 


43.9 




4 q+0.4 


372ti 


5.0+0.3 


767t|* 


17.5+0.8 


10 


PG 2112+059 


2007 May 21 


85.7 


87.5 


66.3 


494+26 


5.8+0.3 


456t^5 


5.2+0.3 


1289+4? 


19.4+0.6 


10 


PG 2112+059 


2007 Nov 05 


50.2 


49.0 


43.0 


190+17 


3.8+0.3 


200+17 


4 1+0-4 
^•-■--0.3 


5901^0 


13.7+0.7 


10 



^ The exposure times and photon counts are obtained after screening the data for flaring. 

^ For details about the source and background extraction regions sec ^ The errors on the source counts were computed by propagating the 
asymmetric errors on the total and background counts using the approach of iBarlowj | |2004D . The total and background count errors were obtained 

from Tables 1 and 2 of[Ge hrcls ( 198^. 

r eferences: (1) ITmanishi Terashimal (1200411 : f 2l ISchartel et all J2005I) ; (3 ) iMiller et a l.' ('2006°); (4 ) IGrupe et al.l lf200^ ; f5') lPage et al.l ||2011D : 
(•6) |Braito et all j2004D : (Tl lRuiz et~aI1 j2007l) ; fS^ IShemmer et al.| ( |2005l) : O^ lBallo eraTj (2011): ao) ISchartel et al.| ( f2010|) . 
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however, in cases where these regions were covering other 
X-ray sources, the background spectra were extracted 
from source- free circular regions with radii of 150" in 
the vicinity of the source. The numbers of counts in 
the source regions at energies of 0.5-2 keV are presented 
in Table [2j The RMFs used for spectral fitting were 
obtained from the public calibration database. The 
ARFs were generated with the pcarf routine; this rou- 
tine does not include aperture corrections, and therefore 
the ROSAT fluxes were aperture-corrected using the in- 
tegrated ROSAT PSPC PSF. 

2.4. ASCA observations 

For the ASCA observations of our sample we used xs- 
ELECT to extract the spectra from the previously pro- 
cessed clean events files publicly available in the Tar- 
tarus (Version 3.1) database. We have extracted spec- 
tra from only the SIS cameras. This was based on the fact 
that the SIS cameras have a wider range of energy sensi- 
tivity and better angular resolution than the GIS propor- 
tional counters. For the SIS, the counts from each quasar 
were extracted from circular regions with aperture radii 
of 3.2'. The background spectra for the SIS observations 
were extracted from source-free polygons covering most 
of the nominal chip at distance > 4.2' from the source 
center. The background subtracted SIS count rates in 
the 0.6—9.5 keV band are presented in Table 2] 

2.5. BeppoSAX and Einstein observations 

For the BeppoSAX and Einstein observations of our 
sample we collected the count rates from already pub- 
lished works. The count rates for the BeppoSAX ME GS 
observations were obtained from 'Braito et al.J (I2004J) in 
the case of Mrk 231 and from Brandt ct al. iflQQQf ) in 
the case of GSO 755. The published count rates were 
extracted in the 1.8-10 keV band using circular regions 
with radii 2' and 3' for Mrk 231 and GSO 755, respec- 
tively. The count rates for the Eins tein IPG observation s 
of our sources were obtained from IWilkes et all ()1994D . 
These count rates were extracted in the 0.16-3.5 keV 
band using circular regions with 3' radii. We used the 
Portable, Interactive M ulti-Mission Simulator software 
(PIMMS; lMukailll993f ) to derive X-ray fluxes from the 
published BeppoSAX and Einstein count rates. 

2.6. New X-ray detection of IRAS 14026+4341 

In this work we report for the first time an X-ray de- 
tection of IRAS 14026+4341, a highly polarized BAL 
quasar that has bee n stu died in detail at other wave - 
lengths fe.g.. Boroson fc Meverslll992l:lHines et al.ll2001[ ). 
We detect this source with a statistical significance of 
« 4(7 in the M0S2 observation performed on 2005 Nov 
26, ft! 3cr in the MOSl observation performed on 2006 
Jan 16, and ~ 5(t in our new AGIS-S observation per- 
formed on 2010 Jul 28. The significances of the detec- 
tions of IRAS 14026-1-4341 were calculated by estimating 
the Poisson probability of obtaining at least the observed 
number of counts based on the expected background. 
This probability provides a statistical test of the null 
hypothesis that the observed photon counts solely orig- 
inate from the background. For the M0S2 observation 

ftp://legacy.gsfc.nasa.gov 
http;//tartarus. gsfc.nasa.gov 



performed on 2005 Nov 26 there are 21 counts detected 
with a predicted background number of 6.63, and there- 
fore the null probability is 6.6 x 10^^; i.e., we reject the 
null hypothesis with a significance of ~ 4.4tT. For the 
MOSl observation performed on 2006 Jan 16 there are 
19 counts detected with a predicted number of 7.45, and 
therefore the null probability is 2.8 x 10""* (significance of 
« 3.4(t). For the Chandra observation there are 6 counts 
detected with a predicted background number of 0.22, 
and therefore the null probability is 1.3 x 10^^ (signifi- 
cance of « 5.1ct). The detections of IRAS 14026+4341 
in the MOS cameras are close to the detection threshold, 
and therefore different choices of energy range and/or 
source and background extraction regions cou ld result in 
a non-detection of the source as in the case of lRuiz et al.l 
(2007^ . For the observations of IRAS 14026+4341 we 
have used our default annular background region. As an 
additional test of the detections of IRAS 14026+4341 in 
the MOS cameras, we select two 80" radius source-free 
circular regions in the vicinity of the source as the back- 
ground. Using these background regions for the M0S2 
observation performed on 2005 Nov 26, the predicted 
numbers of background counts are 7.98 and 9.66. Thus 
the detection significances are ~ 3.7(T and « 3.1(7. For 
the MOSl observation performed on 2006 Jan 16 the 
predicted numbers of background counts are 6.41 and 
7.76; therefore the detection significances are « 3.9(7 and 
« 3.3(7. The detection of IRAS 14026+4341 in its Chan- 
dra observation has « 5(7 significance independent of the 
selection of the background region. For example, using 
as background regions two 20" radius source-free circular 
regions in the vicinity of the source, the expected num- 
bers of background counts are 0.24 and 0.32. Thus, the 
detection significances are « 5.1(7 and rs 4.7(7. 

3. SPECTRAL ANALYSIS 

The spectral analysis of our sample was undertaken 
with the main goal of searching for long-term flux vari- 
ability within the set of X-ray observations of each 
source. Our objective is not to perform a detailed spec- 
tral modeling of each source in our sample but to ob- 
tain robust basic flux measurements. Spectral analysis 
was performed using xspec v. 12.0; the fltting bands used 
were 0.5-8 keV for the Chandra observations, 0.5-10 keV 
for the XMM-Newton MOS observations, 0.3-10 keV for 
the XMM-Newton pn observations, 0.6-9.5 keV for the 
ASCA SIS observations, and 0.5-2 keV for the ROSAT 
PSPG observations. The spectral flts were performed 
in the fu fl fitting ba nd of each instrument using the C- 
statistic (|Gashlll9^79| ) . The C-statistic may not be appro- 
priate to utilize when performing background subtrac- 
tion. To use the C-statistic properly, we fit the back- 
ground spectra with a flat response using the CPLINEAR 
model. The background model is scaled and subtracted 
when we fit the source spectra. 

For Chandra we performed fits using each of the new 
and archival observations. For XMM-Newton we per- 
formed fits using each of the EPIG cameras with detected 
spectra separately; we also performed joint fits using the 
three EPIG cameras together. The fits from each sepa- 
rate EPIG camera are consistent with the joint fits, and 

This model is included in the ACIS-EXTRACT package 
Proos et al.ll20T0h . 
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TABLE 4 

Log of ASCA Archival Observations 











EXP. 


time(ks)^ 


COUNT RATE(x10'^) 




OBJECT NAME 


OBS 


DATE 


SISO 


SISl 


SISO 


SISl 


ref.'' 


IRAS 07598+6508 


1996 


Oct 


29 


37.2 


37.2 


<2.3 


<3.3 


2 


PG 0946+301 


1998 


Nov 


12 


70.1 


64.8 


1.3+0.5 


<1.0 


3 


PG 1001+054 


1999 


Nov 


12 


38.0 


40.1 


<1.4 


<1.3 


1 


Mrk 231 


1994 


Dec 


05 


19.4 


18.2 


9.9+1.2 


7.4+1.2 


2 


Mrk 231 


1999 


Nov 


10 


86.4 


90.0 


11.2+0.6 


9.1+0.5 


4 


PG 1700+518 


1998 


Mar 


24 


18.3 


14.7 


<2.1 


<2.1 


2 


PG 2112+059 


1999 


Oct 


30 


31.1 


28.4 


20.2+1.2 


13.3+1.0 


5 



Note. — The SIS count rates are for the 0.6—9.5 keV band extracted from a circular re- 
gion of radius 3.2'. Count rate upper limits are at the 3cr significance level. The count rates 
and exposure times for all the observations with exception of the one of IRAS 07598+6508 
have been obtained selecting source and background regions from the cleaned event files of 
the NASA Tartarus database (version 3.1; http://astro.ic.ac.uk/Research/Tartarus). The up- 
per l imits and expos u re tim es for the observation of IRAS 07598-|-6508 have been obtained 
from lGallagher et all 119991) . 

^ The exposure times and photon counts are obtained after screening the data. 

^ references: (1) th i s wor k; (2) IGallagher et al.| 119991) : (3) iMathur et al.| (|2000D ; 

(4) IMalonev ReynoTdl ( f2000|) (5) IGallagher et al.| (f200T|) . 
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Fig. 1. — Number of sources versus rest-frame time span of X-ray coverage (left panel) and number of X-ray observations (right panel). 



therefore we present only the joint fits in this work. The 
same conclusion was obtained for the ASCA SIS cam- 
eras, and therefore we present only the joint fits of the 
SISO and SISl cameras in this work. 

For deriving fluxes and luminosities of each source we 
used a simple power-law model including intrinsic ab- 
sorption in the cases where significant absorption was 
found (> 99% level of significanc e); in all cases we 
includ ed Galactic absorption from iDickev fc Lockm"anl 
(1X9901). Galactic absor ption column densities (from 
IDickev fc LockmanI [l990l ) were obtained at the optical 
coordinates of the sources (see Tables 1 and 2) through 
the use of the HE AS ARC Nu tool." The power-law 
photon index (F) was free to vary for those Chan- 
dra, XMM-Newton, and ASCA observations where the 
signal-to- noise ratio (S/N) of the observation was greater 
than seven. The value of S/N is estimated from the total 

http;/ /heasarc. gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh. pi 



number of counts divided by the average of the upper 
and lower error bars; the total photon counts with their 
error bars are found in Tables [H [3] and [H In the cases 
where S/N < 7 and for the ROSAT observations,^^ we 
fixed the photon index to the weighted average if there 
are other observations of the source with S/N > 7 or 
otherwise to 1.8. It is important to emphasize that we 
are not correcting for intrinsic absorption when we de- 
rive fluxes. This allows us to make reliable estimations 
of the absorbed fluxes and luminosities regardless of the 
model that we are using since these quantities are rather 
insensitive to spectral model choice. Given the low S/N 
of most of the spectra analyzed in this work, a simple 
or absorbed power-law gives an acceptable fit for almost 
every observation; however, in a few cases of high S/N 

1* The 0.5-2 keV band observations with ROSAT cannot be 
used to give a reliable estimation of the full-band (0.5-8 keV) pho- 
ton index. 
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TABLE 5 

Best-Fit Power-Law X-ray Spectral Parameters 

















OBJECT NAME 


OBS. 


DATE 


(10^2 cm-2) 


r 


C-stat(DOF) 


Instrument 


IRAS 07598+6508 


1991 


Mar 17 


<0.15 


2.71 


22.0(22) 


PSPC 




2000 


Mar 21 


<0.52 


2.71 


7.0(8) 


ACIS-S 




2001 


Oct 24 


<0.06 


2.71+0.19 


301.3(345) 


EPIC 




2010 


Jun 18 


<0.12 


2.71 


36.4(27) 


ACIS-S 


PG 0946-1-301 


1998 


Nov 12 




1.48 


303.7(321) 


SIS 




zUiU 


Jan 11 


"■3-5.3 


1 /iQ_l_n Qfi 






PG 1001-1-054 


1992 


May 17 


<0.19 


2.04 


11.5(9) 


PSPC 




2003 


May 04 


<0.01 


2.04+0.15 


557.0(325) 


EPIC 




2010 


Jan 11 


<0.27 


2.04 


4.3(11) 


ACIS-S 


PG 1004-f 130 


2003 


May 04 


<0.01 


1.39+0.04 


1204.4(1339) 


EPIC 




2005 


Jan 05 


<0.15 


1.03+0.06 


521.8(545) 


ACIS-S 


Q 1246-057 


1991 


Dec 21 


<8.1 


1.58 


22.3(28) 


PSPC 




1993 


Jan 05 


<3.6 


1.58 


15.5(12) 


PSPC 




1993 


Jun 22 


<2.6 


1.58 


39.0(66) 


PSPC 




2000 


Feb 08 


o 0+3.5 
d.»_2 g 


1.58 




ACIS-S 




2001 


Jul 11 




1.83+0.20 


441.8(584) 


EPIC 




2004 


Jul 15 


<1.2 


1.40+0.17 


542.9(606) 


EPIC 


Mrk 231 


1991 


Jun 7 


<0.01 


1.03 


121.4(103) 


PSPC 




1994 


Dec 05 


<0.10 


1.00+0.22 


251.4(304) 


SIS 




1999 


Nov 10 


<0.02 


0.94+0.10 


814.8(882) 


SIS 




2000 


Oct 19 


<0.00 


1.03+0.06 


731.8(629) 


ACIS-S 




2001 


Jun 07 


<0.00 


1.06+0.04 


1919.7(1663) 


EPIC 




2003 


Feb 03 


<0.00 


1.07+0.06 


719.9(603) 


ACIS-S 




2003 


Feb 11 


<0.00 


1.05+0.06 


847.9(582) 


ACIS-S 




2003 


Feb 20 


<0.00 


0.95+0.07 


694.0(605) 


ACIS-S 




2010 


Jul 11 


<0.03 


0.89+0.20 


189.2(211) 


ACIS-S 


IRAS 14026+4341 


2005 


Nov 26 


<4.5 


1.80 


24.2(18) 


EPIC 




2006 


Jan 16 


<1.5 


1.80 


30.3(18) 


EPIC 




2010 


Jul 28 


<113 


1.80 


8.6(5) 


ACIS-S 


CSO 755 


2001 


Jul 30 


<3.4 


J^-OU_Q 32 


47.8(63) 


EPIC 




2001 


Dec 08 


1.20+0.30 


1.83+0.07 


1036.9(1161) 


EPIC 




2001 


Dec 13 


-LI Q 

2.8ti;| 


2.02+0.29 


213.1(238) 


EPIC 




2010 


Mar 21 


<5.3 


1.97+0.46 


79.0(82) 


ACIS-S 


SBS 1542+541 


1993 


Ago 14 


<4.5 


1.74 


13.5(24) 


PSPC 




2000 


Mar 22 


<4.0 


1.55+0.49 


56.6(65) 


ACIS-S 




2002 


Feb 03 


<1.6 


1.73+0.26 


214.7(254) 


EPIC 




2002 


Feb 06 




1.76+0.13 


576.9(647) 


EPIC 


PG 1700+518 


2009 


Dec 19 


<0.04 


1.77 


49.2(44) 


EPIC 




2009 


Dec 31 


<0.08 


1.77 


141.2(124) 


EPIC 




2010 


Jan 02 


<0.12 


-I 77+0.33 
'-0.30 


191.0(195) 


EPIC 




2010 


Jun 26 


<0.11 


1.77 


19.0(13) 


ACIS-S 



PG 2112+059 



1991 Nov 30 


<0.20 


0.91 


39 


.9(46) 


PSPC 


1999 Oct 30 


1.14+0.50 


1.91+0.26 


506, 


.1(599) 


SIS 


2002 Sep 01 


<0.03 


0.84+0.09 


390, 


.8(399) 


ACIS-S 


2003 May 14 


<0.02 


1.59+0.06 


896, 


.1(958) 


EPIC 


2005 Nov 20 


<0.00 


0.84+0.08 


1095 


.1(1023) 


EPIC 


2007 May 03 


<0.02 


0.66+0.11 


583, 


.6(585) 


EPIC 


2007 May 19 


<0.01 


0.69+0.07 


1053, 


.0(1164) 


EPIC 


2007 May 21 


<0.01 


0.62+0.06 


1495, 


.4(1509) 


EPIC 


2007 Nov 05 


<0.02 


0.58+0.09 


779 


.6(891) 


EPIC 



Note. — The best-fit column densities were obtained from a power-law model with intrinsic 
absorption at the redshift of the source. The same model was used to estimate the best-fit photon 
index and C-statistic for the observations of PG 0946-f-301, Q 1246-057, SBS 1542+541, CSO 755, 
and PG 2112-1-059. The best-fit photon index and C-statistic for all the other observations were 
obtained from an unabsorbed power- law model. In all the fits w e have assumed Galactic absorption 
with column densities obtained from fP^ickey &: Lockmanl l |1990D . Errors and upper limits represent 
90% confid ence intervals for each value, taking one parameter to be of interest (AC— 2.71; e.g., 
I Avni|[T97l) . 

^ Intrinsic neutral hydrogen column density in units of 10^^ cm^'^. 
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Fig. 2. — Logarithm of the observed flux in the 0.5—2 keV band (in units of 10~^^ erg cm~^ s~^) versus modified Julian date for the BAL 
quasars in our sample. The circles, squares, inverted triangles, and triangles indicate observations performed with Chandra, XMM-Newton, 
ASCA, and ROSAT, respectively. Given that Mrk 231 did not present variability in the set of Chandra observations performed in 2003, 
for visualization purposes we have plotted these observations as a single point with their mean flux and a typical error bar. The same 
procedure has been applied to the XMM-Newton observations of IRAS 14026+4341, the XMM-Newton observations of CSO 755 performed 
in 2001 Dec, the XMM-Newton observations of SBS 1542+541, the XMM-Newton observations of PG 1700+518, and the XMM-Newton 
observations of PG 2112+059 performed in 2007 May. The upper limits at early times found in the panels of PG 1001+054, PG 1004+130, 
Q 1246—057, and Mrk 231 correspond to the Einstein observations. The shaded areas in each panel correspond to the expected range of 
X-ray fluxes obtained assuming a typical ctox (see Table 5 of Stcffcn ct al. 2006) and a power-law spectrum with F = 1.9. The sources 
marked with * present evidence of variability in at least one of the energy bands considered. The horizontal dotted lines correspond to the 
best-fit constant flux. 
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Fig. 3. — Logarithm of the observed flux in the 2—8 keV band (in units of 10^^* erg cm~^ versus modified Julian date for the BAL 

quasars in our sample. The circles, squares, stars and inverted triangles indicate observations performed with Chandra, XMM-Newton, 
BeppoSAX , and ASCA, respectively. See the Figure [2] caption for additional details about the presentation. 



(S/N > 25) where we find evidence of complexity in tfie 
spectra, we also estimate fluxes and luminosities using 
more complex models (see §4.2p . 

The best-fitted parameters using a power-law model 
are presented in Table [5j There are five sources 
that presented evidence of absorption in their spec- 
tra; these sources are PG 0946+301, Q 1246-057, 
CSO 755, SBS 1542-1-541, and PG 2112-1-059 (see Ta- 
ble El. PG 0946-f301 shows in both Chandra and 
ASCA observations moderate-to-strong absorption with 
column densities ^ 10^^ cm~^. This source presents 
moderate and low S/N in its Chandra and ASCA ob- 
servations, respectively. Therefore, new Chandra and/or 
XMM-Newton observations with exposures > 20 ks are 
needed in order to better constrain the absorption found 
in PG 09464-301. Q 1246-057 shows moderate intrin- 
sic absorption (A'h ~ 10^^ cm^^) as we confirm in 
the Chandra observation performed in 2000 and the 



XMM- Newton observation performed in 2001 (first stud- 
ied bv iGrupe et al.l [20031 ) . CSO 755 presents moder- 
ate intrinsic absorption (TVh ~ 10^^ cm~^) as we con- 
firm in the XMM-Ne wton observations (first studied by 
IShemmer etldl [20051 ) . SBS 1542-^541 shows low intrin- 
sic absorption (TVh ^ 10^^ cm~^) as we confirm in the 
last XMM- Newton observatio n performed in 2002 (first 
studied bv iGrupe etHTI I2003D . PG 2112-1-059 presents 
evidence of moderate absorption only in its ASCA ob- 
servation performed in October 1999 (first studied by 
Gallagher et al. 2001); however, the fiat measured power 
laws in other observations suggest absorption is likely 
present. 

The derived X-ray fluxes in the 0.5-2 keV, 2-8 keV, 
and 0.5-8 keV bands from simple or absorbed power-law 
fits are presented in Table |6| and Figures [21 [3l and [H 
There we have also presented fiuxes from observations 
with BeppoSAX and Einstein; these fiuxes were ob- 
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BAL quasars in our sample. The circles, squares, and inverted triangles indicate observations performed with Chandra, XMM-Newton, 
and ASCA, respectively. See the Figure [2] caption for additional details about the presentation. 



tained from published count rates using the PIMMS tool 
(see We assumed power-law spectra with the av- 

erage parameters (F and A^h) taken from spectral fits 
of the high S/N observations of each source. The er- 
ror bars for the BeppoSAX MEGS observations were ob- 
tained by propagating the uncertainties on the instru- 
m ental background^^ a nd so urce count r ates presented 
in iBrandt et all (jl999| ) and iBraito et ah (2004). The 
Einstein observations of our sources provided only up- 
per limits. We confirmed that the estimates of Einstein 
fluxes with PIMMS give a reliable re s ult by comparing 
our values with those of I Wilkes et al.l ()1994[ ): in general 
the differences are less than 10%. Similarly, our Bep- 
poSAX fluxe s are c onsis tent within 10% w ith those of 
IBrandt et all ()1999f ) and lBraito et al.l (|2004l ). 

Note that we have only estimated fluxes and luminosi- 
ties in the bands where the instruments are sensitive. 
Therefore, in the case of the ROSAT PSPC and Einstein 
IPC observations we only present fluxes in the 0.5-2 keV 
band, and in the case of the BeppoSAX MECS we only 

To estimate the background count rates we used 
the ba ckground files found at the NASA HEASARC site 
( jhttp:/ /heasarc. nasa.gov/docs/sax/archive/calibration_files.htmlJ . 



present fluxes in the 2-8 keV band. 

4. RESULTS AND DISCUSSION 

In this section we describe the main results found in 
the X-ray variability study of our BAL quasar sample. 
In ^4.11 we constrain variability using simple power-law 
models (as described in ^ to estimate fluxes; in ij4.2l we 
constrain variability using more complex models to esti- 
mate fluxes; flnally, in ^4.3! we examine possible physical 
scenarios to explain the results found. 

4.1. Long and short-term variability 

Long-term variability was assessed initially by compar- 
ing the fluxes in the 0.5-2 keV, 2-8 keV, and 0.5-8 keV 
bands for all the detections of each source. In order to 
test the signiflcance of any variability, we calculate for 
each source the statistic where the data points are 
the fluxes of each observation (with their respective la 
errors) and the model is the best-fit constant flux. The 
value provides a statistical test of the null hypothesis 
that the flux of each epoch is equal to the best-fit flux 
(average flux weighted by errors) of all the observations. 

In Table |7| we list the best-fit constant fiuxes for 
each source in each band. In this table we also have 
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TABLE 6 

Derived X-ray Fluxes from Power-Law Fits 



OBJECT NAME^ 


OBS 


DATE 


/0.5-2'' 






/2keV*^ 


Instrument 


iHAb U75yoH-boUo 


1991 


Mar i / 


£jr»_i_r» on 
2.dU±U.o9 






QO_l_n on 
Z.!58±U.99 


PoPO 




1996 


Dec 21 


<5.0 


<1.9 


<7.0 


<5.6 


SIS 




2000 


Mar 21 


^•^^-0.66 




3.03_Q gQ 




ACIS-S 




2001 


Oct 24 


3.20±0.29 


1.20±0.11 


4.40±0.39 


3.55±0.32 


EPIC 




2010 


Jun 18 


2.10±0.38 


0.79±0.14 


2.89±0.52 


2.33±0.42 


ACIS-S 


PG 0946+301 


1993 


Nov 14 


<1.2 






<0.8 


PSPC 




1998 


Nov 12 


0.68±0.42 


8 3±3 6 


8.9it3.9 


<0.6 


SIS 




2010 


Jan 11 


1 43+018 

'■■^'^—0.27 


ii.2t^:t 


12.6l^-J 


0.69±0.11 


ACIS-S 


PG 1001+054* 


1979 


Dec 02 


<14 






<24 


TOO 




1992 


May 17 


^ j^-t-fl fi4 
^•^"-0.62 






« *-t-i 1 


PSPC 




1999 


Nov 12 


<2.5 


<2.3 


<4.8 


< 4.2 


SIS 




2003 


May 04 


3.53±0.36 


3.35+0.34 


6.88+0.70 


5.99+0.61 


EPIC 




2010 


Jan 11 


0.93±0.27 


0.88+0.25 


1.81+0.52 


1.58+0.45 


ACIS-S 


PG 1004+130* 


1980 


May 09 


<7.7 






<20 


IPC 




2003 


May 04 


11.76±0.38 


27.45+0.89 


39.2+1.3 


28.27+0.92 


EPIC 




2005 


Tan 

J dill \JO 


1 n Qo-i-n ^9 


41.7+1.2 


OZ.U^ J- .u 


on 79-1-0 88 
zy. 1 z^u.oo 


ACIS-S 




197E 


Jul 20 


\ 1-^ 






<38 


IPC 


1991 


Dec 21 


9 07-1-0 57 






4.4+1.2 


PSPC 




1992 


Jul 17 


<2.2 






<6.2 


PSPC 




1993 


Tan 

J dill \jo 


Q ic-i-n Qi 






1 Cl-LQ K 
-LZ.U^O.U 


PSPC 




1993 


Tun 22 


9 91 -l-n '^fi 






7.5+1.2 


PSPC 




2000 


Feb 08 


-L. 1 O^U.OU 


4.9+1.0 


6.6+1.4 


1.97+0.41 


ACIS-S 




9nni 

^UU J. 


Till 1 1 


1 70-1-0 1 


9 AO-I-O 18 


A 10-1-0 "^l 


98-1-0 48 


EPIC 




zuu^ 


Tnl 1 ^ 


1 "^0-1-0 10 


99-1-0 9"^ 


A c;o_|_n Qc; 

"rl.OZ^U.Ou 


Q4-I-0 "^0 


EPIC 


Mrk 231 


1980 


Dec 23 


<12 






<33 


IPC 




1991 


Jun 07 


1 o^u.yo 






QA fi-h9 
01. u^z .0 


PSPC 




1994 


Tier 


1 -L.oin 1 .u 


/I r: Q_Lf; 

'4:0-y^u.o 


c;7 Q-LV Q 
( .o^ / .0 




^T^ 

010 




1 QQQ 


INUV _LU 




^O-O^z.y 


uy.u^o.u 


0-I-. 1 ^ 1 . y 


010 




2000 


Oct 19 


10.83±0.34 


41.7+1.3 


52.6±1.7 


29.48±0.93 


ACIS-S 




2001 


Jun 07 


12.07it0.35 


44.3+1.3 


56.3±1.6 


32.26±0.94 


EPIC 




2001 


Dec 29 




46.0+2.5 






MECS 




2003 


Feb 03 


10.98±0.33 


39.7+1.2 


50.7±1.6 


29.21±0.89 


ACIS-S 




2003 


Feb 11 


10.88±0.34 


40.7+1.3 


51.6±1.6 


29.29±0.91 


ACIS-S 




2003 


Feb 20 


10.29±0.33 


44.1+1.4 


54.4±1.8 


29.09±0.95 


ACIS-S 




2010 


Jul 11 


10.62±0.96 


49.8+4.5 


60.5±5.5 


40.0±2.8 


ACIS-S 


IRAS 14026-1-4341 


1992 


Jul 18 


<1.0 






<2.1 


PSPC 




2005 


Nov 26 


1.28±0.44 


1.69+0.58 


3.0±1.0 


2.73±0.94 


EPIC 




2006 


Jan 16 


0.84±0.42 


1.11+0.55 


1.95±0.97 


1.79±0.89 


EPIC 




2010 


Jul 28 


1 97+0.64 
-0.48 


1 fl7+0-84 
J--O'_0.64 




2 7+1'* 
^•'-1.0 


ACIS-S 




1999 


Feb 02 




11.2+3.8 






MECS 






Tni 

J Ul ou 


Q 0-1-1 9 




OQ 7-m 7 

zo. / ^0. 1 


Q9 K-l-r. n 


EPIC 




2001 


UcC UO 


Q OQ-l-0 9Q 


1 Q 1 0-1-0 49 


00 1 Q-l-0 71 

ZZ. J.O^U. 1 -L 


31.1+1.0 


EPIC 




2001 


Dec 13 




11.1+1.1 


19.6+1.9 


18.8+1.9 


EPIC 




2010 


Mar 91 
iVldi ^ J. 


c Q4-|-n S7 


8.9+1.1 




on c:_|_o 


ACIS-S 


SBS 1542-1-541 


1993 


Ago 14 


6.2±1.3 






15.4+3.2 


PSPC 




2000 


Mar 22 


4 01 -l-n 62 


9.0+1.4 


1 0+9 


10.3+1.6 


ACIS-S 




2002 


Feb 03 


K ^4-1-0 "^0 




13.9+1.3 


19.6+1.8 


EPIC 




2002 


Feb 06 


c 01 _|_n on 


Q 7Q-I-0 44 


14 fiO-l-O 7^ 


22.1+1.1 


EPIC 


pn 1 7004- "^1 >i 


1 QQ1 


1^ u uy 


^1 9 






^9 ^ 


pcpp 




1998 


Mar 94 


<3.0 


^4 9 


^7 9 

\ ( .z 


<6.4 


SIS 




2009 


Dec 19 


u. oozcu. -L y 


1 1 f^-un 9R 




1 77-1-0 40 


EPIC 




2009 


Dec 31 




n fiC-l-O 90 
u.uo^u.zu 


1 1 Q-l-O 'i4 


1 OR-l-0 '^0 


EPIC 




2010 


Jan 02 


0.76±0.16 


1.04±0.21 


1.80±0.37 


1.61+0.33 


EPIC 




2010 


Jun 26 


0.66±0.18 


0.90±0.25 


1.56±0.43 


1.39+0.38 


ACIS-S 


PG 2112+059* 


1991 


Nov 30 


7.0±1.0 






19.6+2.9 


PSPC 




1999 


Oct 30 


21.2±3.8 


62±11 


83±15 


77+14 


SIS 




2002 


Sep 01 


3.17±0.14 


15.81±0.69 


18.98±0.83 


8.94+0.39 


ACIS-S 




2003 


May 14 


17.35±0.74 


30.8±1.3 


48.1±2.1 


42.6+1.8 


EPIC 




2005 


Nov 20 


1.25±0.07 


6.26±0.34 


7.51±0.40 


3.54+0.19 


EPIC 




2007 


May 03 


1.88±0.13 


11.98±0.86 


13.86±0.99 


5.42+0.39 


EPIC 




2007 


May 19 


1.67±0.08 


10.26±0.47 


11.93±0.55 


4.81+0.22 


EPIC 




2007 


May 21 


1.74±0.06 


11.86±0.44 


13.60±0.51 


5.04+0.19 


EPIC 




2007 


Nov 05 


1.21±0.07 


8.73±0.52 


9.94±0.59 


3.53+0.21 


EPIC 



Note. — Fluxes are obtained from the best-fitting parameters for a simple power-law model (most cases) or 
an intrinsic absorbed power-law model (for observations as indicated in the previous ialjle). In both cases, fixed 
Galactic absorption was also included. The fluxes are not corrected for intrinsic absorption. 
^ The sources marked with * present evidence of variability in at least one of the energy bands considered. 
^ Observed fluxes in the 0.5—2, 2 — 8, and 0.5—8 keV bands in units of lO'"^^ erg cm~^ s~^. 
*^ Observed flux densities at rest-frame 2 keV in units of 10~^^ erg cm~^ s~^Hz~^. 
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TABLE 7 

Tests for Variability: Constant Flux Fits 





Fluxes [10 ^'^ erg cm 


-2 g-l]a 




X^/dof 




% 


SIGNIFICANCE 


OBJECT NAME 


0.5—2 keV 


2-8 keV 


0.5-8 keV 


0.5-2 keV 


2-8 keV 


0.5-8 keV 


0.5-2 keV 


2—8 keV 


0.5—8 keV 








POWER-LAW 


FITS 










IRAS 07598+6508 


2.79 


1.05 


3.84 


5.6/3 


5.7/2 


5.6/2 


86.8 


94.3 


93.9 


PG 0946+301 


1.21 


10.46 


11.63 


2.3/1 


0.5/1 


0.7/1 


86.9 


50.8 


58.1 


PG 1001+054 


1.93 


1.89 


3.89 


33.4/2 


32.8/1 


32.7/1 


>99.9 


>99.9 


>99.9 


PG 1004+130 


11.27 


32.50 


44.74 


3.0/1 


87.7/1 


44.2/1 


91.5 


>99.9 


>99.9 


Q 1246-057 


1.54 


2.76 


4.38 


15.7/5 


12.5/2 


4.2/2 


99.2 


99.8 


87.5 


Mrk 231 


11.05 


42.56 


53.45 


18.1/8 


20.0/8 


13.1/7 


98.0 


99.0 


93.0 


IRAS 14026+4341 


1.14 


1.50 


2.63 


0.6/2 


0.6/2 


0.6/2 


26.4 


26.4 


26.4 


CSO 755 


8.82 


12.54 


21.42 


5.7/3 


15.2/4 


9.6/3 


87.2 


99.6 


97.8 


SBS 1542+541 


5.51 


8.75 


14.32 


6.8/3 


0.1/2 


0.7/2 


92.4 


4.3 


28.1 


PG 1700+518 


0.68 


0.93 


1.61 


2.5/3 


2.5/3 


2.5/3 


53.1 


52.8 


52.9 


PG 2112+059 


0.00 


0.00 


0.00 


0.0/0 


0.0/0 


0.0/0 


>99.9 


>99.9 


>99.9 






POWER-LAW FITS WITH 


CROSS-CALIBRATION 


UNCERTAINTIES 






IRAS 07598+6508 


2.63 


0.98 


3.62 


3.4/3 


3.5/2 


3.4/2 


66.5 


82.5 


81.9 


PG 0946+301 


1.18 


10.35 


11.47 


2.1/1 


0.4/1 


0.6/1 


84.9 


48.6 


55.6 


PG 1001+054 


1.61 


1.56 


3.20 


20.0/2 


19.9/1 


19.9/1 


>99.9 


>99.9 


>99.9 


PG 1004+130 


11.29 


31.82 


44.05 


0.3/1 


7.4/1 


3.8/1 


38.9 


99.4 


94.9 


Q 1246-057 


1.59 


2.81 


4.45 


9.5/5 


6.8/2 


2.8/2 


91.1 


96.7 


75.8 


Mrk 231 


11.11 


43.58 


54.39 


2.5/8 


3.5/8 


2.2/7 


3.9 


9.9 


5.2 


IRAS 14026+4341 


1.13 


1.49 


2.62 


0.6/2 


0.6/2 


0.6/2 


25.3 


25.3 


25.3 


CSO 755 


8.24 


11.52 


19.92 


2.1/3 


7.0/4 


4.2/3 


45.7 


86.4 


76.4 


SBS 1542+541 


5.24 


8.74 


14.04 


3.9/3 


0.1/2 


0.3/2 


72.9 


2.2 


14.4 


PG 1700+518 


0.67 


0.92 


1.60 


2.2/3 


2.1/3 


2.2/3 


46.2 


45.9 


46.0 


PG 2112+059 


1.59 


9.58 


11.18 


146.0/8 


99.5/7 


107.1/7 


>99.9 


>99.9 


>99.9 




POWER-LAW FITS WITH CROSS-CALIBRATION 


UNCERTAINTIES AND 


UPPER LIMITS 




IRAS 07598+6508 


2.62 


0.98 


3.61 


3.4/4 


3.5/3 


3.4/3 


50.7 


67.8 


66.9 


PG 0946+301 


1.00 






3.8/2 


85.0 






PG 1001+054 


1.59 


1.51 


3.12 


21.8/4 


20.2/2 


20.1/2 


>99.9 


>99.9 


>99.9 


PG 1004+130 


10.45 






9.2/2 


99.0 






Q 1246-057 


1.58 






11.5/7 






88.4 






Mrk 231 


11.04 






4.4/9 






12.0 






IRAS 14026+4341 


0.86 






3.0/3 






61.4 






PG 1700+518 


0.68 


0.94 


1.62 


3.0/5 


3.0/4 


3.0/4 


30.5 


44.1 


44.0 




FITS 


USING COMPLEX MODELS 


WITH CROSS-CALIBRATION UNCERTAINTIES 






Mrk 231 


10.63 


48.40 


59.13 


4.0/8 


1.2/7 


1.1/7 


14.7 


0.9 


0.8 


CSO 755 


8.60 


11.14 


19.80 


1.0/3 


3.1/3 


2.0/3 


20.1 


63.4 


42.8 


PG 2112+059 


1.39 


10.53 


11.91 


108.1/8 


72.6/7 


81.6/7 


>99.9 


>99.9 


>99.9 



^ Best-fitting constant flux (weighted average) from the Chandra, XMM-Newton, BeppoSAX , ASCA, and ROSAT detections of each source. 
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TABLE 8 

Average S/N, Expected Fraction of Variable Sources, and Observed 
Minimum-to-Maximum Flux Ratios. 



OBJECT NAME^ 


(S/N) 






[ /mi„ J 0.5-2 




L/minJo.5- 


IRAS 07598+6508 


6.6 


50 


0.20 


1.5+0.4 


1.5+0.4 


1.5+0.4 


PG 0946+301 


4.7 


27 


0.20 


2.1+1.4 


1.4+0.7 


1.4+0.7 


PG 1001+054* 


7.8 


69 


0.20 


3.8+1.3 


3.8+1.3 


3.8+1.3 


PG 1004+130* 


45.1 


2079 


0.80 


1.1+0.2 


1.5+0.2 


1.3+0.2 


Q 1246-057 


8.8 


86 


0.20 


2.4+0.8 


2.0+0.5 


1.6+0.4 


Mrk 231 


29.8 


917 


0.43 


1.2+0.2 


1.3+0.2 


1.2+0.2 


IRAS 14026+4341 


2.1 


7 


0.02 


1.5+1.0 


1.5+1.0 


1.5+1.0 


CSO 755 


17.2 


313 


0.43 


1.3+0.3 


1.8+0.4 


1.5+0.4 


SBS 1542+541 


13.4 


193 


0.43 


1.6+0.5 


1.1+0.2 


1.1+0.2 


PG 1700+518 


4.8 


28 


0.20 


1.7+0.7 


1.7+0.7 


1.7+0.7 


PG 2112+059* 


25.2 


660 


0.43 


17.5+4.1 


9.9+2.3 


11.1+2.6 



^ The sources marked with * present evidence of variabiUty in at least one of the energy bands 
considered. 

^ Expected number of counts in the Chandra full band in order to achieve the average S/N of the 
set of observations of each source. 

^ Variability detection probabilities obtained from the expected mean number of counts in the 
Chandra full band. These probabilities are obtained from column 3 of this table and from Tables 1 
and 3 of lGibson &: BrandtJ 120121) . 



TABLE 9 

Derived X-ray Fluxes from Complex Models Fits 



OBJECT NAME* 



OBS. DATE 



fo.i 



/0.5-8 



2kcV 



Instrument 



Mrk 231 



1991 Jun 07 
1994 May 12 

1999 Nov 10 

2000 Oct 19 

2001 Jun 07 
2003 Fob 03 
2003 Feb 11 
2003 Feb 20 
2010 Jul 11 



11.0+1.6 
10.8+1.7 
13.8+1.5 
10.4+1.1 
11.5+0.5 
10.22+0.66 
10.13+0.60 
9.60+0.66 

1-2.1 
-1.5 



10.3^ 



46.0+7.2 
49.1+5.4 
48.9+5.2 
51.9+2.3 
44.1+2.8 
48.8+2.9 
49.6+3.4 
HO 

-7 



49^ 



56.8+9.0 
62.9+6.9 
59.4+6.4 
63.4+2.8 
54.3+3.5 
58.9+3.5 
59.2+4.1 

-12. 
-9 



59"* 



14.8+3.9 
34.7+5.5 
23.7+2.6 
20.5+2.2 

17.88+0.79 
18.3+1.2 

14.93+0.89 
15.2+1.1 

23.5^^5 



PSPC 
SIS 
SIS 
ACIS-S 
EPIC 
ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S 



CSO 755 



PG 2112+059* 



2001 Jul 30 


8.5+1.7 


12.5+2.6 


21.0+4.3 


26.7+5.5 


EPIC 


2001 Dec 08 


9.14+0.43 


12.46+0.59 


21.6+1.0 


27.4+1.3 


EPIC 


2001 Dec 13 


8.6+1.2 


10.4+1.4 


19.0+2.6 


17.4+2.4 


EPIC 


2010 Mar 21 


7.1+1.6 


8.3+1.9 


15.4+3.5 


22.3+5.1 


ACIS-S 


1991 Nov 30 










PSPC 


1999 Oct 30 


31.3+4.2 


63.8+8.6 


95+13 


75+10 


SIS 


2002 Sep 01 


3.06+0.32 


14.7+1.5 


17.7+1.8 


7.14+0.74 


ACIS-S 


2003 May 14 


16.3+1.3 


34.7+2.7 


51.0+3.9 


41.2+3.2 


EPIC 


2005 Nov 20 


1.13+0.08 


7.54+0.56 


8.67+0.64 


2.20+0.16 


EPIC 


2007 May 03 


1.62+0.22 


12.6+1.7 


14.2+1.9 


3.67+0.50 


EPIC 


2007 May 19 


1.54+0.18 


10.7+1.2 


12.3+1.4 


4.06+0.46 


EPIC 


2007 May 21 


1.44+0.11 


12.38+0.91 


13.8+1.0 


3.62+0.26 


EPIC 


2007 Nov 05 


1.12+0.13 


9.1+1.1 


10.2+1.2 


3.02+0.36 


EPIC 



Note. — Fluxes are obtained from the best-fitting parameters of the ionizcd-absorbed power-law model 
for CSO 755 and the thermal plasma + two absorbed power-laws model in the case of Mrk 231 and 
PG 2112+059. The fluxes are not corrected for intrinsic absorption. 

* The sources marked with * present evidence of variability in at least one of the energy bands considered. 
^ Observed fluxes in the 0.5 — 2, 2 — 8, and 0.5 — 8 keV bands in units of lO^^'^ erg cm~^ s^^. 
^ Observed flux densities at rest-frame 2 keV in units of 10""^^ erg cm~^ s~^Hz~"^. 



listed the values of for our model to test vari- 
ability. We find that six out of eleven sources show 
evidence for variability (significance > 99%; see Ta- 
ble [7]). The sources with potential variability are 
PG 1001+054, PG 1004-f-130, Q 1246-057, Mrk 231, 
CSO 755, and PG 2112+059. Notice that for most of 
these sources the apparent variability is found between 
observations performed with different X-ray instruments 
(e.g., see Figures [21 [3l andH]). The only exception is 
PG 2112-1-059 where variability is independently found 
in the XMM-Newton observations (see Figures [2 [3l S] 



and Table Therefore, in order to assess the vari- 
ability reliably we must take into account uncertainties 
associated with cross calibration between the different 
X-ray missions used in this work. From studies of differ- 
ent astrophysical sources that can be considered "stan- 
dard candles" in the X-ray regime the main conclusion is 
that, in general, the cross-calibration errors in the X-ray 
fluxes between Chandra, XMM-Newton, BeppoSAX , 
ASCA, ROSAT, and Einstein s hould b e < 10% (e.g 
Snowdenll200"2l:lKirsch et al.ll2005l: Ncvalai nen et al.ll2010l: 



Weisskopf et al.ll2010HTsuiimoto et al.ii20n[) . Therefore. 
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as a second test of the variability found we add in quadra- 
ture to the error bars of each measured flux a 10% error. 
The new constant flux estimates with their values and 
significances are presented in Table [71 We find that three 
(PG 1001+054, PG 1004+130, and PG 2112+059) out 
of the five variable candidates show significant variability 
(at a level > 99%) in the measured fluxes. 

We have made a basic comparison of the fraction of 
X-ray variable BAL quasar s in our sample with the l arger 
Chandra quasar sample of iGibson fc Brandtl ()2012l) . To 
make an appropriate basic comparison with this sam- 
ple, we obtain the average S/N of the observations of 
each source in our sample (see Table [S]). These S/N val- 
ues are converted to an equivalent number of counts in 
the Chandra full band (see Table [8]). We assume that 
the error s in the Chand ra counts are given by Tables 1 
and 2 of iGehreli (|1986f ). We proceed by obtaining the 
expected variability detection probability for each source 
based on the cou nt-dependent fract i on of variable sources 
found in §3.2 of IGibson fc Brandtl (l20ll). For this pur- 
pose, we have divided the sources in IGibson fc Brandtl 
(|2012[) into four count ranges: 1-10, 10-100, 100-1000, 
and 1000-10000 counts. The fractions (variability detec- 
tion probabilities) of detectably variable sources in each 
range are 0.02, 0.20, 0.43, and 0.80. For each source 
we select one of these four values based on its equiva- 
lent average counts in the Chandra full band (see Ta- 
ble [5]) . Adding the probabilities found in Table [51 the 
expected number of variable sources is « 3.5. Therefore, 
the fraction of X-ray variable sources appears compara- 
ble to_tliat_ex£ected for no rmal quasars. As is shown 
in IGibson fc Brandt (|2012[ ). the fraction of detectably 
variable sources in X-rays has a first-order dependency 
on the S/N of the observations of each source. This 
is because the sensitivity to measure any change in the 
flux increases with S/N. Additionally, the UV luminosi- 
ties of our so urces are broadly compar able to those of 
the sample of IGibson fc Brandtl (|2012D . although both 
sets of sources span large ranges of luminosity (see Fig- 
ure [SJ. The mean luminosity of our sample ((L2500) = 
30.95 +0.31) is similar to the m ean luminosity of the sam- 
ple of IGibson fc Brandtl (I20TI) ((L2500) = 30.55 ± 0.04). 
Therefore, any secondary luminosity dep endency (e.g., 
iNandra et all il997t [Manners et "all 120021) ^'^ of the vari- 
ability should not fundamental l y affe ct our basic com- 
parisons with IGibson fc Brandtl ([2012D . 

The ratio of maximum-to-minimum fluxes of the 
variability is 3.8 ± 1.3, 1.5 ± 0.2, and 9.9 ± 2.3 for 
PG 1001+054, PG 1004+130, and PG 2112+059, respec- 
tively. These ratios have been calculated in the 2-8 keV 
band since in this band all the sources present variability. 
Additionally, the shortest rest-frame timescales of signif- 
icant variability were approximately 5.8, 1.4, and 0.5 yr 
for PG 1001+054, PG 1004+130, and PG 2112+059, re- 
spectively. Since in general there is not extensive cov- 
erage of sources with a single X-ray mission (with the 
exception of PG 2112+059 and Mrk 231), the flux uncer- 
tainties of our variability analysis are greater than 10%. 
Hence, for most of the sources we are not sensitive to 
amplitude changes < 1.4.^^ In general the amplitude 

Not cl ear relation was fo u nd be tween luminosity and X-ray 
variability in IGibson fc Brandti ([2013). 

^"^ Assuming two flux measurements with a fractional error of 
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Fig. 5. — Logarithm of the monochromatic luminosity at 
rest-frame 2500 A versus redshift. Monochromatic luminosities 
used here have units of erg s~ ^Hz~^. Black-flUed circles are data 
from IGibson & Brandt I I2012I) . Our sample is represented with 
open triangles. 

threshold for variability detection is higher than 1.4, es- 
pecially for observations with less than ~ 100 counts, 
for which errors in the flux are above 10%. As found 
in Table [5] the maximum-to- minimum flux amplitudes of 
the sources that do not present statistically significant 
variability are 10-100% in the 2-8 keV band. 

Notice that in the variability searching above we have 
not taken into account flux upper limits. In general, as 
seen from Figures [H [31 and [H the typically loose flux 
upper limits do not present any new useful information 
regarding variability; specifically, the upper limits are 
usually above or not far below the best-fit constant fiux. 
As one basic illustrative test to see if the upper limits 
give any new information about variability, we assume 
that the upper limits are detections with fluxes given by 
one half of the threshold and error bars that are three 
tenths of the upper- limit value. Under this assumption 
we test for variability by calculating the fits of constant 
flux in Table [71 The test with the upper limits presented 
in Table [71 is with 10% cross-calibration errors taken into 
account. From this test we find that the upper limits are 
possibly giving new information about variability in just 
one case. This case corresponds to the Einstein upper 
limit for PG 1004+130. From Figure H and Table [3 it 
is clear that the upper limit increases the significance 
of the soft-band variability of PG 1004+130. Since this 
source has already shown variability in the hard band 
and possibly full band, adding the upper limit does not 
contribute any material change to the conclusions of the 
variability analysis. 

We also searched for rapid variability in the 
high signal-to-noise ratio (S/N > 25) Chandra, 
XMM-Newton, ASCA, and ROSAT observations by ap- 
plying the Kolmogorov-Smirnov test to the photon ar- 
rival times. No significant variability (at a level of sig- 

10% each, our variability test requires a maximum-to-minimum 
flux ratio of > 1.4 in order to get a > 99% level of significance. 

If we assume that the flux below the upper-limit has a uniform 
distribution, the mean value is half the threshold and a is ^1/12 Ri 
0.28 of the upper-limit. 
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nificance >99%) in the full band of each instrument was 
found. 

As illustrated with the shaded regions in Figures [U [31 
and[4l the BAL quasars in our sample are usually X-ray 
weak by factors of « 2-100 relative to expectations for 
non-BAL quasars (see §4.31 for further discussion). Thus, 
any substantial removal of the shielding gas responsible 
for this X-ray weakness (see 21) from our line of sight 
should lead to large X-ray flux changes — changes that are 
larger than the upper limits that we typically place above 
upon allowed X-ray variability. The fact that such large 
X-ray flux changes appear rare over 3-30 yr suggests that 
the shielding gas is fairly stable, even on timescales suf- 
ficiently long to allow significant physical rearrangement 
and accretion-disk rotation. This fact should be consid- 
ered in future m odeling of BAL q uasar X-ray absorption 
variability (e.g.. ISirn et al.ll20100 as well a s in interpre- 
tatio n of UV BAL variability results (e.g . . iGibson et al. 
2008L [20101 : iCanelhiPO et~al1[20Tll [20ll iFiliz AkeTaL 

mm . 

4.2. Flux variability in complex spectra 

For low S/N observations it is not possible to probe 
effectively spectral complexity that deviates from the 
power-law models used in ^ However, in cases where 
the S/N is sufficiently high we can quantify these devi- 
ations. We address this issue by grouping with a min- 
imum of ten counts per bin the high S/N observations 
(S/N > 25) and then performing F-tests to assess the 
improvement of the fits from a power-law through the 
use of more complex alternative models. 

The alternative models selected are used to probe com- 
plex absorption and to corr oborate previous studies of 
sources in our sample (e.g., 'Malonev & Revnolds' '2000'; 
Gallagher et al. 2002, 2004; Schartcl et al. 2005, 2010; 
Miller et al.ll200ai .^^ Through this method we find that 
in three sources at least one of the alternative models 
improves the fit with significance > 99% in every high 
S/N observation; these sources are Mrk 231, CSO 755, 
and PC 2112+059. In the case of Mrk 231 there were 
two models (models 6 and 7 in Footnote [191) that pro- 
vide improvement of the fits in every high S/N obser- 
vation. These models are thermal plasma -|- two ab- 
sorbed power-laws and thermal plasma -I- reflection spec- 
trum. Note that when we use a model with two absorbed 
power-laws, each power-law component has the photon 
index fixed to F = 2 and the column density free to vary. 
The power-law that is heavily absorbed represents the 
intrinsic X-ray emission and the second power-law with 
less absorption repre s ents scattered em ission (see, e.g., 
[Gallagher et al.ll2002t [Braito et ai][2004[ ). In the case of 
CSO 755 an ionized absorbed power-law model improves 
the f it in the only high S /N observation (in agreement 
with [Shemmer et al.[ [20051 ). In the case of PG 2112-^059 
there were two models that show improvement of the fits 

^® These models are (1) absorbed power-law with iron emis- 
sion line, (2) ionized-absorbed power-law (xSPEC model ABSOai), 
(3) partially covered absorbed power-law, (4) reflection spec- 
trum (xsPEC model PEXRAv), (5) thermal plasma (xSPEC model 
vraymond) -|- power-law, (6) thermal plasma -|- two absorbed 
power- laws (scattered spectra), (7) thermal plasma -|- reflection 
spectrum, (8) reflection spectrum from ionized media (xSPEC model 
PEXRiv), and (9) power-law with absorption from a blueshifted iron 
line. 



in every high S/N observation. These models are ther- 
mal plasma -I- two absorbed power-laws and reflection 
spectrum from ionized media. 

Based on these fitting results, we have recalculated 
the fluxes in each observation and band for Mrk 231 
and PG 2112-1-059 using the thermal plasma -I- two ab- 
sorbed power-law model, and for CSO 755 using an 
ionized-absorbed power-law model. The fluxes ob- 
tained using these new models are presented in Table [9l 
and plotted in Figure IHl The models used to obtain 
the fluxes in Figure |6l (and Table [91 involve more pa- 
rameters than a simple power-law; therefore, the error 
bars on the fluxes in this figure are in general some- 
what larger than those obtained in Figures [H [31 and SI 
The flux-variability analysis for Mrk 231, CSO 755, and 
PG 2112+059 using complex models is consistent with 
that using a simple power-law model (see Table [3 and 
^4.ip . Additionally, the weighted average fluxes obtained 
for Mrk 231, CSO 755, and PG 2112+059 through the 
use of more complex models are within 20% of those ob- 
tained using a power-law model (see Table[71). Therefore, 
as expected, our main variability results appear robust 
to alternative modeling of the spectra. 

4.3. Spectral variability 

As expected our sample of BAL quasars is generally 
X-ray weak when compared with RQQs of similar UV 
luminosities (see Table [TOl and Figures [2l [31 HI and ^ . 
On average our sample has luminosities at rest-frame 
2 keV that are ~ 20 times weaker than expe ctations 
for RQQs ([Steffen et al.l[200l lJust et al.l[2007l) . With 
the exceptions of CSO 755 and SBS 1542+541, which 
are as X-ray luminous as typical RQQs, all of our sam- 
ple is > 3(7 below the aox — ^uv a nd L2keV ~ ^ uv re- 
lations (see Table [101 and Table 5 of lSteflFen et all [2006'). 
In Figure [3 we have plotted ctox and L2kGV versus Luv 
for our sample. For the cases of X-ray variable sources 
we have marked the states of maximum and minimum 
X-ray flux in the aox — -^uv and L2koV — ^uv diagrams 
assuming that the UV luminosity is constant. The as- 
sumption that the UV luminosity is constant is required 
given that the X-ray observations do not have simulta- 
neous optical/UV follow-up in general. However, in gen- 
eral it is expected that the luminosity in the UV will 
vary much less than that in the X-rays. For example, 
the XMM-Newton observations of PG 2112+059 present 
maximum-to-minimum flux ratios in the X-ray s of sa 6 
and in the optical of « 1.1 (jSchartel et al.|[2bl0f) . 

To perform a basic test for X-ray spectral variabil- 
ity, we utilize the effective photon index from a simple 
power-law fit, recognizing that this photon index may not 
have a direct physical interpretation in cases of complex 
spectra. Using the same variability test as in §4.11 we 
conclude that changes in F have significances > 99% for 
PG 1004+130 and PG 2112+059.^1 For PG 2112+059 

^'^ The error in the estimated fluxes using the thermal plasma + 
two power-law model depends mainly on the normalization of the 
two power-laws (npLi and npL2), and therefore we estimate the 
error on the flux in a given band (/bp) by A/bp = | ^np^i l'^"PLi + 
hf^|AnpL2. 

S"PL2 ' 

We have taken into account the cross-calibration errors in V 
which are estimated to be < 5% fe.g.. ISnowden|[2002l : iKlrsch et al.l 
[20051: iTsuTimoto et al.|[20T^) . 
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Fig. 6. — Logarithm of the observed flux (in units of 10~^* erg cm~'^ versus modified Julian date for Mrk 231, CSO 755, and 

PG 2112+059. The fluxes were obtained using alternative models to a simple power-law. The three upper, middle, and lower panels 
correspond to fluxes in the soft (0.5-2 keV), hard (2-8 keV) and full (0.5-8 keV) bands, respectively. The circles, squares, inverted 
triangles, and triangles indicate observations performed with Chandra, XMM-Newton, ASCA, and ROSAT, respectively. See the Figure[2] 
caption for additional details about the presentation. 



the changes in T are confirmed with significant variabil- 
ity (> 99% significance) in the complex fits hard-to-soft 
flux ratios (/2-8//0.5-2) (see, e.g., Table[9]). In the case 
of PG 1001-^054, the low S/N of the Chandra observa- 
tion does not allow us to constrain the photon index well, 
and therefore we cannot tell if there is any change in the 
shape of the spectra. We conclude that the flux variabil- 
ity of PG 1004-^130 and PG 2112+059 shows associated 
spectral variability. 

From our spectral analyses alone it is not possible to 
explain the extreme weakness of the X-ray luminosi- 
ties of our sources; however, this is likely associated 
with complex absorption. As already noted in ^ al- 
most half of our sample (flve out of eleven) presents di- 
rect evidence of absorption. Additionally, the flat effec- 
tive spectral indices measured for most of the sample 
could be an indication of absorption. For example, as 



described in iMiller et al.l (|2006[ ). the hard spectrum of 
PG 1004+130 could imply partially covering absorption. 
Furthermore, based on our spectral analyses it is likely 
that the observations of Mrk 231 and the weak states 
of PG 2112+059 correspond to epochs of heavy absorp- 
tion (A^H > 10^'^ cm~^). Therefore, the complexity in 
the observed spectra of these sources could be the re- 
sult of heavily absorbed emission plus scattered emission 
(see m.'2\i . In the case of Mrk 231 this picture has been 
shown to be plausible in the joint spectral analysis of 
XMM-Newton and BeppoSAX obs ervations at e nergie s 
between 0.2-50 keV performed bv iBraito et al.l ()2004D . 
The spectral softening of PG 2112+059 with increasing 
flux might be explained by an increase in the observable 
X-ray emission due to a decrease in the absorption (see 
^i^. Alternative interpretations of the physical picture 
for PG 2112+059 (e.g., a reflection dominated spectrum) 
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Fig. 7. — Logarithm of the monochromatic luminosity at rest-framo 2500 A versus Qox (upper panel) and logarithm of the monochromatic 
luminosity at rest-framo 2 keV (lower panel). Monochromatic luminosities used here have units of erg Hz~^. Black-filled circles are 
data from Stcffcn ct al. (200(|); arrows indicate X-ray upper limits. Dashed lines are the best-fit linear relations for their combined sample. 
Our sample is represented with filled and open triangles for sources with and without variability, respectively. For variable sources we have 
plotted the state of minimum and maximum flux joined with a dotted line. PG 1004-1-130 does not show variability in soft X-rays, and 
therefore its 2 ke^V luminosity is presented as a single point. 

ulation on rest-frame timescales of 3-30 yr. For seven out 
of the eleven sources we have also obtained and reported 
on ne^w Chandra observations suitable for detecting any 
strong X-ray variability. Our main conclusions are the 
follo'wing: 



can be found in lSchartel et all ()2010[ ). The ASCA obser- 
vation of PG 2112-1-059 presents a typical quasar spectral 
shape ■with absorption that is moderate enough so the 
majority of the X-ray emission can get through the ab- 
sorber. Therefore, the X-ray emission of PG 2112-1-059 
in its brightest state is likely to be dominated by un- 
absorbed emission over scattered emission. The harden- 
ing of the spectra -with increasing flux for PG 1004-1-130 
is not easily explained ■with absorption alone; its origin 
could be more com plex due to X-ray jet-linked emission 
(for details sec Mill er et al.l[2006.) . Future X-ray observa- 
tions, especially ■with high sensitivity above 10 keV such 
as those with NuSTAR and ASTRO-H, will help to clar- 
ify the physical picture of these sources. 

5. SUMMARY AND FUTURE WORK 

We have analyzed Chandra, XMM-Newton, Bep- 
poSAX, ASCA, ROSAT, and Einstein observations of 
eleven BAL quasars to perform the first systematic in- 
vestigation of the long-term X-ray variability of this pop- 



1. We find that three (PG 1001-^054, PG 1004-1-130, 
and PG 2112-1-059) of the eleven sources sho^w 
significant flux variability (at a confidence level 
> 99%). The maximum observed amplitude of 
this variability in the 2-8 keV band is a factor of 
3.8 ± 1.3, 1.5 ± 0.2, and 9.9 ±2.3 for PG lOOH-054, 
PG 1004-1-130, and PG 2112+059, respectively. 
These three sources sho'w detectable variability on 
rest-frame timescales do'wn to 5.8, 1.4, and 0.5 yr, 
respectively. See §4.1 and §4.2. 

2. The flux variabflity of PG 1004+130 and 
PG 2112+059 has associated spectral variability 
(the data for PG 1001+054 do not allo-w effective 
investigation of spectral variability) . This spectral 
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TABLE 10 

Optical and X-ray Properties of Surveyed BAL Quasars 



OBJECT NAME 


Log L2500'' 


Log L2-10'' 


Log L2kov'' 


ALog L2kov'* 


Log /2_8'= 


ALog f2-s^ 




AOox'' 


ixxi-vo Uioyo-roouo 




.0 1 


9A 1 S 
Z^. ±0 


— Z.ZO 








— U.OO 


PG 0946+301 


31.68 


44.74 


25.42 


-1.86 


-12.99 


-0.23 


-2.40+0.08 


-0.66 


PG 1001+054 


29.69 


42.34 


24.22 


-1.65 


-13.81 


-1.70 


-2.10+0.07 


-0.64 


PG 1004+130 


30.58 


43.82 


25.61 


-0.89 


-12.50 


-0.62 


-1.91+0.05 


-0.32 


Q 1246-057 


32.10 


44.94 


26.65 


-0.93 


-13.55 


-0.46 


-2.09+0.05 


-0.29 


Mrk 231 


29.04 


42.36 


24.08 


-1.33 


-12.36 


-1.02 


-1.90+0.04 


-0.53 


IRAS 14026+4341 


30.50 


42.78 


24.80 


-1.65 


-13.83 


-1.61 


-2.19+0.10 


-0.62 


CSO 755 


32.32 


45.88 


27.63 


-0.11 


-12.94 


0.25 


-1.80+0.05 


0.02 


SBS 1542+541 


31.94 


45.51 


27.32 


-0.15 


-13.06 


0.21 


-1.77+0.05 


-0.00 


PG 1700+518 


30.77 


42.48 


24.48 


-2.16 


-14.03 


-2.11 


-2.42+0.06 


-0.80 


PG 2112+059 


31.29 


43.99 


25.40 


-1.61 


-13.02 


-0.99 


-2.26+0.04 


-0.57 



^ Logarithm of the monochromatic luminosity at rest-frame 2500 A (with units erg s~^ Hz~^). These values were computed from the flux 
densities at rest-frame wavelength 2500 A using the corrected for galactic reddening flux densities of the closest optical magnitude and 
extrapolating a power- law with a ^ -0.5 fe.g.. |Vanden Berk et al.||2001<) . 

^ Logarithm of the 2 — 10 keV band luminosity (with units erg s~^). For each source this quantity was obtained from the weighted- aver age 
luminosities of the Chandra, XMM-Newton, and ASCA observations in the 2 — 10 keV rest-frame band. The luminosities have been corrected 
for Galactic absorption. 

^ Logarithm of the monochromatic luminosity at rest-frame 2 keV (with units erg s~^ Hz""^). For each source this quantity was obtained 
from the weighted- average flux of the Chandra, XMM-Newton, ASCA, and ROSAT observations at 2 keV rest-frame. 

The difference between the measured and the predicted log L2kcV for RQQs, based on the lJust et al relation: 
log L2kcv - 0.709 X log L2500 + 4.822 . 

^ Logarithm of the flux in the 2 — 8 keV band (in units of lO"^'^ erg cm~^ For each source this flux was obtained from the weighted- aver age 

flux of the Chandra, XMM-Newton, BeppoSAX , and ASCA observations in the 2—8 keV band. 

^ The difference between the measured and the predicted flux in the 2—8 keV band for RQQs. The predicted fluxes in the 2—8 keV band arc 
obtained from the expected fluxes for RQQs at 2 keV rest-frame Uust et al.ll2007f) assuming a power-law spectrum with T — 1.9. 

S Ttie optical-to- X-ray power-law slope a^x = log(/2kev//2500A)/log(i^2kev/'^2500A)- 

^ The difference between the measured aox value and the predicted aox for RQQ, based on the lJust et al.l (120071) relation: 

aax = -0.140 X : 



i2500A+ 2.705. 

variability may be induced by absorption changes, 
although the hardening of PG 1004+130 when it 
brightens is puzzling in this regard. See §4.3. 

The eight sources in our sample without signif- 
icantly detected X-ray flux variability are con- 
strained by our data to have varied by less than 
10-100%. See §4.1. 

We do not find that BAL quasars exhibit excep- 
tional long-term X-ray variability when comparing 
the measured variability for our whole sample to 
that of non-BAL quasars. Strong changes in the 
X-ray absorbing shielding gas, owing to physical re- 
arrangement or accretion-disk rotation, appear rel- 
atively rare (although some changes are present). 
Since the shielding gas is critical for setting the ion- 
ization level of the wind that produces UV BALs, 
our results have implications for modeling BAL 
variability. See §4.1. 

We report for the first time an X-ray detection of 
the highly polarized and well-studied BAL quasar 
IRAS 14026+4341. We detect this source with a 
statistical significance of ~ 5a in our new Chandra 
observation. See §2.6. 



Chandra and XMM-Newton observations. Improving the 
number of detected counts per epoch or extending the 
the bandpass coverage to higher energies will enable bet- 
ter physical understanding when absorption changes are 
detected; missions including NuSTAR, ASTRO-H , and 
Athena are required for such observations. 
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This first systematic investigation has significantly im- 
proved constraints upon the long-term X-ray variability 
of BAL quasars in general. However, further improve- 
ments are needed in the number of objects monitored, 
the number of monitoring epochs, the number of detected 
counts per epoch, and the bandpass coverage. Improv- 
ing the numbers of objects and epochs will set better 
constraints upon the frequency of significant absorption 
changes; this could be achieved with, e.g., additional 
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